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Abstract 

We compute the pressure of the standard model at high temperatures in the symmetric phase 
to three loops, or to 0{g^) in all coupling constants. We find that the terms of the perturbative 
expansion in the SU(2) + Higgs sector decrease monotonically with increasing order, but the large 
values of the strong coupling constant Qs and the Yukawa coupling of the top quark gy make the 
expansion in the full theory converge more slowly. The final result is observed to be about 10% 
smaller than the ideal gas pressure commonly used in cosmological calculations. 
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1 Introduction 



The most fundamental quantity for the thermodynamics of a finite temperature gauge theory is its free 
energy. Its partial derivatives give measurable values of thermodynamic quantities and its singularities 
indicate possible phase transitions. Extensive work has been devoted to the study of the free energy 
(F = fV = —pV) of QCD and electroweak theory in various temperature ranges and with various 
levels of rigor. 

For QCD and low T, (T < T,, ^ Aqcd) one can use hadron gas models, the region around T = Tc 
can and has to be studied with first principle lattice Monte Carlo methods [1,2] and at large T the 
method with controllable accuracy is the effective theory one [3,4]. This is based on asymptotic 
freedom and on separating the relevant scales ttT, the electric scale toe = gT and the magnetic scale 
uiM — g^T. The computation of the coefficients of the expansion in g has a long history: the results of 
orders g'^ [5], g^ [6], g'^lng [7], 5^ [8], g^ [4,9] and g^\ng [10,11] are known. The calculation has also 
been generalized to the case when the chemical potentials associated with quarks are nonzero [12]. 
This expansion has several interesting finite-T-effccts: odd powers of g appear, the logarithmic terms 
are logs of the ratios of the matching scales, log(T/mE) and log{mE/mM), and, finally, the coefficient 
of the g^ term is not perturbatively calculable [13], since vacuum diagrams of all orders in the loop 
expansion contribute to it. This is due to the fact that the effective theory of the magnetic sector is 
confining. However, the coefficient of the g^ term can be determined by a combination of numerical 
and several involved analytic computations [14]. One then has the pressure as a well defined expansion 
in the coupling constant g determined in the MS scheme. 

While both QCD and the electroweak sector of the standard model are gauge field theories, there is 
a crucial difference between them: the Landau pole of QCD corresponds to a length scale 1 / Aqcd ~ 
10~^^ m, the length scale of a nucleon, while for the electroweak theory 1/Aew ~ 10^ m, comparable 
to the radius of the Earth. Thus, while confinement effects are important in the QCD case; near Tc, 
rendering perturbative calculations unreliable, they are negligible in the electroweak case and therefore 
it is, at least in principle, possible to apply perturbative methods down to Tc and even below that. 
Due to this there has been an extensive amount of work devoted to calculating the properties of 
the electroweak phase transition using, for example, perturbative l-loop [15-18] and 2-loop [19-21] 
effective potential calculations. Those methods are reliable only for small Higgs masses and the 
c;omplete solution of the problem required first a perturbative matching of the full 4-dimensional 
theory to an effective 3-dimensional theory [22] and then numerically solving the phase diagram from 
the effective theory using lattice Monte Carlo techniques [23-27] . The phase diagram was observed to 
have a first order line which ends in a 2nd order critical point of Ising universality class [28] . Similar 
techniques have been used to solve the phase diagram also when the external U(l) magnetic field [29] 
or the chemical potentials related to the baryon and lepton numbers [30] are nonzero. The phase 
diagram has also been solved with numerical studies of the full 4-dimensional theory [31]. Grand 
unified theories have been studied in [32,33]. 

The structure of the phase diagram and the properties of the phase transition (latent heat, interface 
tension, correlation lengths, order parameter discontinuities) only depend on the discontinuities of the 
free energy, not on its value. In the effective theory approach there is thus an important theoretical 
step missing, the matching of the value of the free energy itself. For QCD this problem was formulated 
and solved in [4, 10] . The purpose of this paper is to do the same for the standard model. At the same 
time we obtain the pressure to order g^. 

To this end we have to carry out a number of rather extensive computations. First, we compute the 
3-loop free energy in the fuU 4d standard model by evaluating the 3-loop finite-T sum-integrals in the 
MS scheme. Schematically, we need F = l+g"^ + g^{l/ e+1) , where the 5^ coefficient contains 1/e terms 
due to infrared (IR) divergences which then cancel against the ultraviolet (UV) divergences of the 
effective theory. The UV divergences of the full theory are cancelled by the standard renormalization 
procedure. Second, we determine the 2-loop screening masses of the gauge fields Aq, Bq in the full 4d 
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theory in the MS scheme. Schematically, we need = g^{l + e) + g^- We then repeat the same for 
the fundamental scalar mass, which is present already on the tree-level but gets thermal corrections. 
Unlike the screening masses, it has additional divergences, so we need = —i''^+g'^{l+e)+g'^{l/e+l). 
Finally, we compute the 3-loop free energy in the 3d effective theory. Schematically, we need / = 
+ g^nn? + gi,m. 

Our final result will have some qualitative differences to the pressure of hot QCD. There is a mass 
scale V independent of T in the Lagrangian which gives rise to terms such as v^T'^, not present in 
PQCD • Another difference is related to renormalization of the fundamental scalar mass in 3d effective 
theory. In a gauge theory with no scalar particles all the parameters of the 3d theory are finite to 
order g'^, whereas in this theory the mass of the fundamental scalar contains divergences requiring 
renormalization. We also have terms of order g^ In g, which cancel in QCD. In addition, the different 
mass scales tod, fn'^ and lead to terms of type g^\nmi/m2, where are some combinations of 
these masses. 

The paper is organized as follows: in section |21 we explicitly define the theory we are working with 
and fix various conventions. In section we briefly review the method of dimensional reduction. 
Sections ^ El and El contain the calculations. Finally, in section [3 we discuss the result, the details of 
which are given in Appendices fXI - IdI 



2 The basic setting 

The theory we consider is the SU(3)c x SU(2)i x U(l)y standard model with np = 3 families of 
fermions and ns = 1 fundamental scalar doublets, and the quantity we will evaluate is the pressure 
of this theory at high temperatures. The theory is specified by the Euclidean action (in the units 

h = c=l) 

S = '^'^ j ^'''^'^ 

+IlIPIl + eRpCR, + qLpqL + urPur + dRpdR + igy {qLT^^*tR ~ ifl(«'*)V2gi) , (2) 



where 



G% = d^A^^d^Al + ge-'-A'^At, F^, - d^B, - d,B^, 
= d^C^-d^C^+gsr'^C'^C', 

pen. = Ifiidf^eR + ig'B^eR) 



PlL = l^id^QL^fA'I^T-qL-^B^qL-igsG-^T-qL), 



2ig' 

pUR = 7p ( df^UR —Bf,UR - igsC^T^UR 

pdR = l^{d^dR+'-^B^dR-igsC;^T-dRY (3) 
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Here A^, 



is the fundamental scalar doublet 



and are gauge bosons of weak-, hyper- and strong interactions, respectively; $ 
II and en arc the left-handed lepton doublets and the right- 
handed lepton singlets (wrt. weak charge), and Ql, ur and da are the left-handed quark doublets 
and the right-handed up and down -type quark singlets. Only the Yukawa coupling for the top 
quark is taken into account. Summation over different families is assumed. Also, d — 3 — 2e in 
dimensional regularization, fi, 1/ = 0,...,d. The gamma matrices are defined in Euclidean space so 
that {7^,7^.} = 2Sfj,„, {75,7/i} = and Tr 757^7^7p7<^ a e^^p^. The color indices are a = 1, c^a 
for the weak interaction and a = 1, — 1 for the strong interaction. The different group theory 
factors for SU(A^) with generators T"^ are defined as: 



FO- = Tr T^T" 

ah ^ace fhce 



CpSij 
dA 



(5''^ dp 



Sii. 



(4) 
(5) 



For SU(2) with = r"/2 they are Tp 1/2, Cp = 3/4, Ca 2, = 3 and rfp = 2. 

The momentum integrations are done using dimensional regularization for both IR and UV di- 
vergences. The dimensionful parameter is chosen according to the MS scheme, which amounts to 
replacing the scale parameter ^ by 



-1/2 



(6) 



All couplings are implicitly scaled to their 4d (e — 0) dimension with /i, so that e.g. 
where g is the coupling in the 4 — 2e-dimensional Lagrangian, [g] = e. We use the Feynman gauge 
— 1) for the gauge particle propagators at all stages of the calculation. The final result should, 
of course, be gauge independent, since pressure is a physical quantity, but we have not checked this 
explicitly. 

The theory contains six couplings that run with the renormalization scale: gauge couplings g', 
g and gs, the fundamental scalar quartic self-coupling A and its mass parameter v^, and gy- The 
counterterms in the standard model can be found in [19]. However, the terms proportional to A were 
neglected there, while the terms including g''^ were dropped in [22], so for completeness we list the 
running of the needed parameters here: 



A(A) 
9sW 



A(m) + 
5'(m) - 



1 

8^ 
1 

8^ 



-3CFff' - 3jg'^ + N,gl + 2{dF + l)x) u'' In -, 
4 J ji 



3 Cp 2 /2 
9 + n-r—r9 9 



2dF + V 



hg'^ + (8 + 2dF)A2 - 3g^ + 2N,Xg^) In -, 



1 

8^ 
1 

8^ 
1 



= 5?(M) + ^ 



9l 



{N^ - 1) g', - y 
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12' 



9y In 



A 



-yCA + ^(A^c + l)^Tp + irpns) / In -, 



-yiV. + -npj,^(M)ln-. 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 



Note that A as we have defined it differs from [19] by a factor 6. Numerically, we fix the values of these 
couplings at the scale /x = mz according to their tree-level relation to different physical parameters: 
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V (mz) 
gyimz) 



'■Hi 



'2 



(mz) = 



mw) , 



V2 

g^(mz) = 4V2G^m 
as{mz) = 0.1187, 



(13) 

(14) 
(15) 



where is the unknown mass of the Higgs boson, mw = 80.43 GeV, mz = 91.19 GeV and 
nit — 174.3 GeV are the masses of the W and Z bosons and the top quark, respectively, and = 
1.664 • 10~^ GeV~^ is the Fermi coupHng constant. We always assume that ng = 1 and np = 3 unless 
stated otherwise, but use the general form to keep track of different contributions. Note, in particular, 
that the result will not be valid for ns > 1, since the mixing of scalars is not taken into account. We 
employ a power counting rule A ~ g'^ ~ ~ ^ 9^ ^^'^ assume the temperature to always be 
so high that the relation v'^ < g^T"^ applies. 

The physical observable we are studying is the pressure, defined by 



PiT) 



T 

lim — In 



VAVxjjV^V^ exp {-S). 



(16) 



It is normalized such that the (real part of the) pressure of the symmetric phase vanishes at T = O.'^ 
The purpose is to calculate the pressure up to, and including, order (7^(l+ln g)T'^, employing the power 
counting rules above. This amounts to calculating 1-, 2- and 3-loop vacuum diagrams contributing to 
the pressure. These, together with their symmetry factors, are listed in Appendix^ Other interesting 
variables, such as entropy and energy densities s(T) and e(T), can then be evaluated using standard 
thermodynamic relations, s(r) = dp/dT, e{T) = Ts(T) - p{T). 



3 Separation of scales 

In this Section we will shortly review the rationale of dimensional reduction applied to our case. A 
more complete treatment can be found, e.g., in [22]. 

A straightforward approach in analytic calculations is to use perturbation theory in evaluating 
the quantities one is interested in. Since we are working in a temperature range where the gauge 
couplings, due to asymptotic freedom, are small, one could naively hope this to be a consistent 
procedure. However, in practice the straightforward expansion in g^ is inhibited by various infrared 
singularities requiring resummations. This in turn leads to the introduction of many different mass 
scales. At high temperature and small coupling the dominant energy scale is the temperature T, while 
the electric and magnetic scales gT and g^T are suppressed by powers of g. The perturbative result 
includes logarithms of all these, making it impossible to choose the UV cutoff in such a way that there 
would not be any large logarithms left. This seems to render perturbation theory unusable. 

The solution, as is well known [3] , is to separate the contributions of different scales into successive 
effective theories, where all the large scales are integrated out one by one. First, we integrate 

p{T) = p^{T) + f PAfcI?AoI?$exp (-5e) , (17) 

where 5'e contains only the static Matsubara modes of the gauge bosons and of the fundamental scalar 
(Higgs) field. The contributions of the nonzero Matsubara modes and fermions to the pressure show 
up as the matching constant pe (Sec^ and in the parameters of 5'e (SecEJ. The spatial (magnetic) 
gauge field components remain massless, while the temporal component gets a thermal mass mn gT. 

^Since the symmetric phase is unstable at T = 0, the pressure there develops an imaginary part when loop corrections 
are calculated. The imaginary part can be related to the decay rate of the unstable phase [34] . 
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The theory defined by Se can then be viewed as a 3d gauge theory with adjoint and a fundamental 
scalar fields. 

The efi^ective theory thus obtained still contains contributions from two scales, gT and g^T, so one 
more reduction step is useful. We integrate out the scale gT and are left with 

p(T) = pE(T)+pM(T) + ^ln y"l?AfcI?$exp(-5M). (18) 

The precise form of the remaining effective theory Sm depends on the conditions of the system. If the 
temperature is much higher than the critical temperature of the system*, then both the adjoint and 
the fundamental scalars can be integrated out since both of them are massive, m ~ gT, and thus the 
remaining effective theory contains only the spatial gauge fields. The only mass scale of the theory 
is then provided by the 3d gauge coupling and is of the order g^T. Consequently, the contribution of 
this theory to the pressure is of the order g^. However, close to the phase transition the fundamental 
scalar, which drives the transition, becomes light and we are not allowed to integrate it out at the 
same time as the adjoint scalars. Then the remaining effective theory contains both the spatial gauge 
fields and the fundamental scalar field, which now has a mass of the order g^l^T . To leading order 
the contribution from this theory to the pressure is then of the order g^l'^T^ — g'^y/g T^. 

In the present paper we will just consider the case when the temperature is much higher than the 
critical temperature, and postpone the study of the case when the system is near the phase transition 
to a later work. The final result of our calculation can then be written as 

piT) = pe{T) + pm{T)+pqcb{T) + 0(/r4), (19) 

where pqcd can be taken from [4,8-10]. One- loop quark diagrams are included in pE, so they must 
be subtracted from pqcd- 

4 Calculation of the pressure p-^ 

The contribution of the nonzero Matsubara modes and fermions to the pressure given by the matching 
constant pb is determined by calculating the path integral in Eq. (|16|l without any resummations 
fADDendix lD.l(l . The calculation involves two diff'erent mass scales, the temperature {2ttT) and the 
mass of the Higgs field (i^^). Since we assume the temperature to always be so high that i/^ < g^T^, 
we can expand the scalar propagator in powers of and keep only terms up to the desired order 
(integration over the scale 27rT is infrared safe and thus pE must be analytic in i/^). The general form 
oi pe(T) can then be written as 

PEiT) = T"^ aEl+ g'^0:EA+ g'^CtEB + ^aEX+ gyCtEY 

+ ^^^^2 (g^OlEAA + g"^a.EBB + {gg' foiEAB + A^OffiAA + Ag^a^AA + W'^OIEBX 
+ gyaEYY + iggY)^aEAY + [g' gY^OlEBY + ^^gyOiEYX 

+ {ggsfaEAs + {g'gsfctEBs + [gYgsfaEY^ 

O^Ei^ + ,, sn ig^OLEAiy + g'^CHEBv + Aa^Ai/ + ffi-aBYi/) 

(47r)^ 

aE.u + T^ ■0{g''). (20) 



(47r)' 



Within perturbation theory there is always a first order phase transition in electroweak theory. 
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The values of all the coefficients oe can be found in Appendix The contribution coming solely 
from QCD is not included. 

All the couplings in the expression above are the renormalized couplings and thus run according to 
Eqs. (|7jl- (|12|l . However, not all the 1/e poles are cancelled by the renormalization procedure as can 
be explicitly seen in the coefficents q;e- The remaining poles correspond to the IR divergences and are 
only cancelled when the contribution from the effective theories to the pressure is taken into account. 

5 Parameters of the 3d theory 

The 3d effective electroweak theory has the general form 

+/i3^^$ASAS + h'^^^^B^B^ - ig3.g^Bo*tASr'^$| , (21) 

where G° = - d,A°; + gge'^'-^^^A^^, = d^Bj - djB^, D,<^> = (9, - ig3T''Af/2 + ig^B,/2)$ and 
D,A^ = diA^ + gae^^M^Ag. AU the couphngs and masses in H21() can be determined to the required 
order in by matching the static Green's functions computed in both the effective and the original 
theory. 

5.1 Coupling constants 

The 2-loop diagrams in the 3d theory are of the order g|m^ ^ 0{g'^). Therefore, the leading order 
results for couplings are enough for our purposes, while the corrections would contribute at g^. At 
tree-level the reduction to the 3d theory only includes scaling the fields by VT, and therefore matching 
the Green's functions gives 

<?i = g'T, g'i = .g'^T, 

A3 = AT, \a = 0(.g*), (22) 
h, = \g^T, h', = y^T. 

The quartic couplings of the adjoint scalars, Xa,b, are not needed at this order. Expressions for 
them can be found in [22]. Note that the relations above hold for (dimensionful) parameters in the 
d-dimensional Lagrangian. The dimensional regularization scale does not need to be the same in 4d 
and 3d, which gives 0{e) corrections to the above matching formulas, e.g. g| = (A/^3)^'^g^T, where A 
and /ia are the 4d and 3d dimensional regularization scales, respectively. At the end of the calculation 
we are going to set A = ^3. When properly renormalized, the 4d theory has only IR divergencies left 
and the scale A should be interpreted as the factorization scale separating the full and the effective 
theory. 

5.2 Mass parameters 

In general, the mass parameters of the effective theory can be found by comparing the poles of static 
propagators in both theories. In the full theory we have for the pole of the propagator 

+ + Tl{k^) + + n(fc2) + Usik^) = 0, (23) 
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at /c^ = k'^ = — m^jf, ko = 0. Here n3(fc'^) is the contribution of n = modes only, and this part is 
also correctly produced by the effective theory, where the same propagator reads 



i^ + n3(fc^) = at k'^ 



(24) 



For m < gT the leading order solution is fc ~ gT. Since n(A:^) has no infrared divergence, we can 
expand it in fc^, 



n(fc2) = n(o) + fc2An(o) 



Up to 0{g'^) Eq. ^ then reads 
k^ 



,2 + n^i^(o) + n'''(o) + n3(fc2) = o. 



and the matching condition can be read from 



d 77(1) 



77(1) 



77(2), 



mi + Usie) = ( 1- (0) ) m^+U' '{0)+n' '{0)+U3ik^) 



(25) 



(26) 



(27) 



givmg 



(28) 



Note that 113 ~ g^mT ~ g^T^ cancels between the two equations, since we only need terms up to g*. 
There are 1/e-divergencies at mlg^ order in the free energy computed in this theory, so we will also 
need the g^e-teims in the masses. 



5.2.1 Adjoint scalar masses 

Applying the procedure described above to the static Aq and Bq propagators gives, after calculating 
all the 2-loop corrections to them. 



(fei + f3E2e + 0{e')) + {pE3 + 0{e)) + 0{g') 



I3exX + 13es91 + 13ey91 + Pe'q'"^ + I3e^- 



(47r)2 



y2 



9" {^'ei + P'E2e + 0{e')) + ^ (/3^3 + 0{e)) + 0{g'') 



(47r) 



/^iJA-^ + f^Esd^ + P'ey9y + + P'eu 



,2 \ 1 



y2 



(29) 



(30) 



The coefRcients (3 are listed in Appendix IbI 

In the above equation g and g' are the renormalized couplings, which run as in Eqs. 1)1 0|l and 
lfTT)l. Substituting these into Eqs. and piljl . we note that at C(e°) all the dependence on the 
dimensional regularization scale A is cancelled by the running of t/'s. 

There are also some electroweak corrections to the adjoint scalar mass of QCD. In addition to pure 
QCD terms already given in [4] we have 



QCD 



2y2y(qcd) 2Cpdp^.g2+2 



11 

'36 



n^g- 



dpgy 



(31) 



and these corrections need to be taken into account in the one-loop term 2m'^^/3TT of the QCD 
pressure. 
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Note that the adjoint scalar masses are finite, unhke the fundamental scalar mass below. This is a 
direct consequence of the fact that the adjoint scalars are actually gauge field components, which have 
no mass renormalizations, and there are no IR divergences in the matching computation. Electroweak 
theory also includes a fundamental scalar field, whose mass contains 1/e poles at two-loop level 
renormalization, since it is not protected by the gauge symmetry. 



5.2.2 Fundamental scalar mass 

In [22] the 3d mass of the fundamental scalar was calculated using effective potential methods, but 
only C'(e'^) terms were given, and g'^ terms were dropped at two-loop level. We calculated this mass 
using the same methods as for the adjoint scalar masses, and included also the 0{g^e) corrections. 

Unlike in the previous section, here we get UV divergencies proportional to that are not cancelled 
by the counterterms of the 4d theory. These are related to the mass renormalization in the 3d theory, 
since the matching procedure gives the bare mass rri^g. It serves as an additional check to calculate 
the 2-loop counterterm dircetly in the 3d theory to see that it precisely cancels the 1/e-terms found 
here. 

For the divergent part we have, substituting directly the correct numerical values for group theory 
factors and setting ns = 1, 



J.2 



81/1 7 ,A 15o,9 
-64^ +64-'^ +32-^"-^ 



4 ^ 



3 \ /2 
4A5 



3A^ 



(32) 



while the counterterm in the 3d theory reads 



1 



(47r)2e 



(33) 



The divergent part is independent of the gauge parameter ^ in covariant gauges. This was expected, 
since the gauge choice in the 3d theory should not depend on that of the 4d theory. 
The finite part gives the renormalized 3d mass. 



1 



9 



(4^) 



9 ■ 



(47r) 



A 



uB 



(An) 



VA 



{Any 



+ T2 [g^i^Ai + f3A2e) + g'^i^Bi + ^32^) + A(/3ai + Px2e) + 5y(/3yi + PY2e) 



9 



AA 



9' 



/4 



2 2 

, 9 9y 

+ , . ^r, PAY 



(47r)' 
9"9'y 



r/3. 



BB 



9^9'\ W « 



A.9 



12 



(47r)2 



(47r)2 



/9. 



All the different coefficients (3xy are given in Appendix IbI 



9I9I r, , Agf- 

^PsY + ,. .^ P\Y 



[Aixf 
(47r)2 



A2 



B\ 



JYY 



(47r)- 



r/3. 



AA 



(34) 



6 Calculation of the pressure py[ 

Computing all the 3-loop vacuum diagrams given by the action (I21|l f Appendix ID.2|) and treating 7713 
as being of order gT produces 



Pm{T) 
T 



An 



dpns {rria + dm 



2\3/2 



16 



3 2m3 I 



1 

An 



-dAm^ 



1 

(4^ 



\—d'p{d'p + l)nsA3m3 — dFdAnsh^m^mu — d-pn^h'^mzm'Y) 
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(47r)3 

+ g'ihmaBBXf + XlmsBxxf + hlm^Bhhf + hlmuBhha + h'^m^iB'f^tif + f^'i'^'u^hht 
+ 535?"i326(m3) + glg'imBb{mB) + glg'im'j^h{m'^) + -^^{dAhsmu + h'^m'j^Y 

47713 

+ dpm^l- + :^— r +33<^ACFdF- In + In 

\Zmu 2mj-| / 6 \mD ma mo / 

+ dpidp + l)Aa(dA^3"iD + /iWd) + glhamuBAha + g'i h'^m'^B'g^^^, + glk'^m'^B^hb 

+ g'.^h3mBBBha+ glh^m^BAhf] ■ (35) 



Constants B^ and the function b{x) are given in Appendix lO Due to divergences in (m| + Smf 



3/2 



'3) 

we have to expand this term in powers of the couphng constants. In higher order terms it is enough 
to use the leading order result 

ml « ml ^ + (^5^ + + \X+ \gl^ (36) 

for the thermal Higgs field mass. 

All the 1/e poles at 0{g^) cancel, and substituting the running parameters of Eqs. ((7t <12|l the g^ 
and orders are seen to be independent of the dimensional regularization scale. The poles at 0{g'^) 
cancel against those in pE coming from the heavy (ttT) modes. 



7 Numerical results 

In this Section we plot the final result given by Eq. I^Hl, into which Eqs. (jSOl, lEHI), (EOl, 

and l|35|) are inserted, for various values of parameters. In particular, we set rriH — 130 GeV, which is 
above the experimental lower limit [35]. Note that vacuum stability considerations lead to a slightly 
higher limit [36], but it turns out that the precise value of rriH does not affect the result much. 

7.1 SU(2) + fundamental Higgs 

Analyzing the result obtained is complicated due to large number of different fields and couplings 
between them in the complete standard model. The total effective number of degrees of freedom of 
the theory is 106.75, but the contribution from the Higgs field that drives the transition to this number 
is just 4 (a complex scalar doublet). Its contribution to the pressure can therefore be expected to be 
small. Also, because the strong coupling constant and the Yukawa coupling of the top quark gy are 
numerically large when compared to the other gauge couplings and to the Higgs self coupling, their 
contribution to the pressure dominates over the contribution coming from the Higgs sector that is 
relevant for the phase transition in the system. It is therefore instructive to consider a simpler SU(2) 
+ Higgs model for which the total number of degrees of freedom is lower and which does not include 
couplings that are not directly related to the phase transition. We achieve this simply by putting 
g'^ = gs — 9y — ^ ^^"i np = in the general result. This theory has also been studied on lattice [37]. 

To lowest order the pressure of this theory is the ideal gas pressure of SU(2) gauge bosons and a 
massless scalar, given by 

Po(T) = ^T4(2dA + 2dF) = ^T^ (37) 
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Figure 1: The pressure of SU(2) + fundamental Higgs theory. The Higgs mass is mn = 130 GeV and 
the W mass is mw = 80 GeV. The critical temperature is Tc — 220 GeV {g^). 

Normalizing the results to po, we plot the pressure of this theory to different orders of the couplings 
in Fig. n The mass of the Higgs boson is taken to be tuh = 130 GeV and the mass of the W boson 
mw = 80 GeV. As can be seen, at high temperatures the introduction of interactions reduces the 
pressure, but since the coupling is small, this effect is small as well. The perturbative expansion is 
well behaved in the sense that the absolute value of each new correction is smaller than that of the 
previous one. This is in contrast to QCD, where the expansion fluctuates around the ideal gas pressure 
unless the temperature is taken to be asymptotically large. 

The result differs from that of QCD also in the sense that there is another manifest mass scale in 
the system, the mass of the Higgs boson. The terms v^T^ and ly'^ in the expansion of the pressure 
become more significant and the pressure deviates from the standard Stefan-Boltzmann law p ^ T'^ 
as the temperature gets smaller. Schematically, the pressure of a gas of massive particles is given to 
leading order by p ~ T'*(l — g^) — T^, where i labels all the particle types in the system and 

the masses are the thermal masses, ~ g^T^ for the temporal component of the gauge bosons and 
J7j2 ^ _jy2 _|_ y2y2 Jqj. ^]-^g Higgs scalar. Thus, as the temperature is lowered, the pressure, normalized 
to the ideal gas pressure of massless particles, behaves as p/po ~ 1 — ff^ + jT"^ . This is seen explicitly 
in Fig. [2 

Although the calculation presented in this paper is not, in principle, valid near the phase transition, 
we have plotted in Fig. ^^Iso the pressure of the broken phase to orders and g^ , corresponding 
to a two-loop calculation. This gives us qualitative understanding about the behavior of the pressure 
near the phase transition. 

The pressure of the broken phase plotted in Fig. ^ is given by being the expectation value 
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Figure 2: The scale dependence of the pressure of SU(2) + fundamental Higgs theory. The temperature 
is fixed to T = 500 GeV, the Higgs mass is nin = 130 GeV and the W mass is mw = 80 GeV. 



of the Higgs field) [19] 
1 



PBp(r, if) 



2 _ . 4 

2 ^ 4 ^ 



9 192^ 24 



^24 ^^"^^^"^ + 3mGs(<y5)^ + ^■mwi'^Y) 



T 

12^ 



+6mvy(<^)3 + 3 ( mw{vf + ^.g'T^ 



3 |^mGs(</3) 

3/2' 



16^ 2 



(38) 



where muiip)^ = 3A(/?^ — i^^, mGs(v)^ = ~ and mi4'(iy9)^ = l/Ag'^ip'^ are the zero temperature 
masses of the particles, and (p = (p{T) is such that dpBP /dip'^ — 0. We can now directly observe that 
there is a temperature where the pressures of the symmetric phase and of the broken symmetry phase 
are equal. Below that the pressure of the broken symmetry phase is bigger and thus the symmetry of 
the theory gets spontanously broken. 

Another interesting question is how the scale dependence of the expansion behaves as higher order 
corrections are added. If the perturbative expansion is well behaved, one expects the scale dependence 
to reduce as more terms are included. This is plotted in Fig. |2] The temperature is fixed to T = 
500 GeV. As can be seen, the result depends very weakly on the chosen renormalization scale: varying 
the scale within A/T ~ 10~^ . . . 10^ changes the result just about two percent. More specifically, 
it is seen that the result up to order 17^ is fairly scale independent. This, however, is a numerical 
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Figure 3: Pressure in different orders of perturbation theory. 

coincidence which stems from the particular values of the parameters , A and v and does not appear 
to have any fundamental reason. The weak scale dependence reappearing in terms of order g'^ and 
seems to support this conclusion. A similar phenomenon was observed in pure gauge theory in [8], 
where the scale dependence of the pressure up to order g^ was seen to be much weaker than expected. 

7.2 The standard model 

Here we plot the pressure of the full theory for realistic values of couplings, using tree-level relations 
between the measured values of toiv, to^, Gf, TUt, ois and the parameters in our result, as shown in 
Eqs. H13l) - (|15|l . The unknown Higgs particle mass is set to the lowest experimentally accepted value 
fUH = 130 GeV. 

For an ideal gas of massless SM particles we would have the familiar Stefan-Boltzmann result 

Po = + 2dA + 2{Nl - 1) + 2dF + 2^nF(dF + 1 + N,{d^ + 2))^ = l06.75^^^ (39) 

which actually is a^i + gluons. This behavior dominates the pressure, so we again divide by po in 
the plots to see the deviations from massless ideal gas. 

The pressure of the full SU(3)xSU(2)xU(l) standard model with 3 families of fermions is plotted 
in Fig. 121 We show the behavior of the result with increasing orders of perturbation theory to find out 
the relative size of corrections. Unlike the SU(2) + Higgs case, the result varies strongly with every 
new order included. This is the known behavior of QCD, and follows from the large values of and 
gy, while the higher order terms in 5, g' and A are small. The 0{g^) correction is still large enough to 
push the line downwards near the phase transition. The relative deviation from the ideal gas pressure 
is of the same order of magnitude as in QCD, which can be explained by the large number of QCD 
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Figure 4: Relative dilFerence between pressure at niH = 130 GeV and mH = 200 GeV. 



degrees of freedom (79 of the total 106.75). As seen above, in the SU(2) + Higgs theory the deviation 
is significantly smaller. We have not plotted the pressure all the way down to T^, where the behavior 
is very singular and the line shoots up to infinity. This stems from the IR divergences in terms like 
TOp/ms, since our assumption m^iT) gT breaks down near and becomes small. 

The effect of varying ran is shown in Fig. 0] where the relative difference between pressure at 
niH — 130 GeV and mn = 200 GeV is plotted. The Higgs particle mass affects the behavior of the 
pressure only very weakly, the only change being a slight and almost constant (times T"*) increase in 
the pressure with increasing tojj. This was expected since the fundamental scalar only represents 4 
of the ~ 100 degrees of freedom. 



8 Conclusions 

In this paper we have calculated the pressure, or the free energy, of the standard model at high 
temperatures to three loops, or to order . The result is about 10% smaller than the ideal gas 
pressure. The effective number of degrees of freedom is thus accordingly reduced from the standard 
value of 106.75 used in cosmological computations. The higher order corrections to the pressure are 
numerically dominated by contributions coming from the strong coupling constant and the Yukawa 
coupling of the top quark. Neglecting them shows that the perturbative behavior of the underlying 
gauge + Higgs theory is good, with the absolute value of each new order in the expansion of the 
pressure being smaller than that of the previous one. This conclusion is supported by the expansion's 
small dependence on the renormalization scale. The large numerical values of gs and gY imply that 
the expansion of the pressure of the complete standard model is not as well behaved. However, even 
then the convergence of the expansion is better than that of the pressure of pure QCD. 

It is possible to extend the expansion presented in this paper by one more order in powers of the 
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coupling constants by using perturbative calculations, to the order g^lng. To evaluate that term 
would require a 4-loop calculation of the vacuum energy densities of the three-dimensional effective 
field theory in Eq. H21|l . However, it is impossible to give an unambiguous meaning for this term until 
the complete term is evaluated as well, a task that already requires the use of nonperturbative 
methods. The convergence of the perturbative expansion is rather fast so one can expect this term to 
be numerically small. 

The calculation performed is valid when the temperature of the system is much higher than the 
critical temperature of the system. However, unlike in QCD, where the coupling constant becomes 
large close to the phase transition and renders the perturbative methods unreliable, in electroweak 
theory it is possible to extend these calculations also to the temperature region close to the phase 
transition. This requires a modified effective 3d theory in which the adjoint scalars Aq and Bq are 
integrated out but the fundamental scalar $ is kept, assuming 7713 is small compared to ttid and mj-,. 
This is left to a forthcoming work [38] . 
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B Matching coefficients 

B.l Coefficients for the adjoint scalair masses 
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Taking into account the different group theoretical factors, the pure gauge and fermionic parts of 
these results agree with [4]. 

B.2 Coefficients for the fundamental scaleir mass 
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D Diagrams contributing to the pressure 

In this appendix we list all the diagrams required for the computation of the pressure. The notation is 
as follows: solid lines represent left-handed fermion doublets and right-handed fermion singlets (thick) 
or just fermion doublets (thin), dashed lines fundamental scalars, dot-dashed lines SU(2) and U(l) 
(thick) or just SU(2) (thin) adjoint scalars, wavy lines SU(2) and U(l) (thick) or just SU(2) (thin) 
gauge bosons and curly lines SU(3) gauge bosons. Dotted lines stand for ghosts. 
Defining the integration measures as 
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The integrals are evaluated in Appendices A and B of [4] and in [39]. 

D.l Diagrams in the full theory 

Here we list the results for the diagrams from the 4d theory. 
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= --CACFdFnFil + N^)g^iD-2) 2A1_2,2 + Mo.o + 2 (£> - 6) J1J1J2 



X hU2,-2 + {D-4 



= ^nUCFTFdF{l + N,fg'' + 



n 2 



4J\A2,-2 + (£> - 4) jVo,o - 4(6 - D)I^l2 




-CACFdFnsg^ 2{D - 2)A^2,-2 



; \-2xi 



D + 2 



Mo,o - 8{D - 2)2f J2 





22 



= CACFdFiD-2)\s9^l!l2 

= i (c^T^d^g'' + ^5'*) ris [4A^2,-2 - A^o.o + 4 (£) - 6) J? J2] 



(158) 



(159) 




4CpTpdp(l + TV,)/ + rfp (^1 + ^ + ^^iVc ) 5'^ 



npns 



o 



-np (2 - DY { -C^dF (1 + A^c) 9^ + 



+lc,dJl + f]g^'-- 



1 17 \ dp A A^c 
- H -/Vc + — H 

2 162 / 32 V 81 



a 



Mi,-i + Mo,o+(li-Ii) Ti 
-\{D- 2f np {N^^ - 1) [c^d^g'gl + ^^g'^gj'^ 



(160) 



(161) 



^'4 



(162) 



A^i,_i + >lo,o + (^Ti - Xi)^ J2 



(163) 
(164) 



-(^D-2)N^[\g-gl^^lg'\l 



7Wi,_i + A1o,o+ (2:1 -T^l2 




{2D - 4) {Nl - 1) 9^5,2 



-1 + Mofi + (^Xi - I2 



/ / 

i i 



2dF {dp + if nsX^Xll2 



(165) 
(166) 

(167) 
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1 x< 



-1 xi 



-OCXj 
< p 

^CX£) 
^d o 



-dpidF + 1) (4CF.g'A + g'^X) nsl^l2 



~24iVeAff^.XiXil2 



= -Z) 



(168) 
(169) 

(170) 

(171) 



icldp/ + ^rfpg'" + JcFdFff^ff'^ ) ns (4Xo,o +Ii'l2) (172) 



0,0 



DN,{?,g^gl+g'^gl)l,IiJ:2 



-Nc {ig^gl + g'hv) (-^0,0 



(173) 



(174) 



(175) 



(176) 



D.2 Diagrams in the effective theory 

The required d = 3 — 2e dimensional integrals have been worked out in the literature except for the 
particular combination of diffenent masses and massless propagators in Eq. (|198(l . The results are 
given here in terms of the integrals /, J, if, A/, iV and L given in the 



2xf « = dpll^ims) 



1 xf } 



1 / V' 



1 / ^/ \ 



-dpidp + l)X3l^{m3) 



-dFdAh3h{m3)Ii{mB) - dFh'^Ii{7n3.)h{'m'Y)) 



(177) 
(178) 
(179) 
(180) 



24 



1 ' \ 



-^^F (Cpgi + ^g'i^ [lUma) + 4:mlMm3)] 



(181) 
(182) 




a 



4 " 



VP 



1 




1 /r-v 

- X 1 • 



CaCf^fSs 



(183) 



-/Vi -1(^3) +ii,-i(?ri3) + Mo,o(m3) - 6/1 (7723) Ji (7723) + -iVo,o(m3) 



+ 6m37Vi,o(m3) - ^m\Mx,Q{mi) + 4m3iVi,i(m3)] 



-iMi,_i - ^M_i,i - i/i Ji - m^Mi,o 



(184) 
(185) 



1113 

-iVi,_i(TOD) + --Li,-i(tod) + -Mo,o(tod) - -/i(tod)Ji(tod) 



5 3 

+ -^0,0(^0) + -mD7Vi,o(mD) - 2mDMi,o(mD) + mD7Vi,i(mD) 



(186) 



1 X 



1 

- X 



1 . 

2"' 




C^rfF - JcaCf^f) 53' + 2CFrfF^53'53' + ^53' 
X [Mi,_i(m3) - 4Mo,o(m3) + 27i(m3) Ji(to3) + %m\Mxfi{mz)\ 



Cld^gi 



(187) 
(188) 

(189) 




o ' \ . I 



C|dF - icACF^F) 53 + 2CFdF^53ff3' + ^Ss' 



dF(dF + l)A^A/o,o(m3) 




3c^ 

= CirfA5'3^^o,o(mD) 



dMo,o(m3) (190) 
(191) 
(192) 



xf V \ 



hldrdAlEims, m^, mn, mn) + h'^dFTEims, ms, rn'o, m^) 
+ 2CFdF^glg3lE{m3,m3,mD,m'c) 



(193) 
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(194) 



C^T^d^gi+'^g'i 



^— I Ji ( Ji + 4m^i<:2) + -jyJ—^ {No,o + SmlNi^o + 16m|7V2,o) 



Aid -I) 



3d — 2 

= -Cld^gi—^ — - [Mo,o + 4m^Mo,i 



(195) 







(196) 



(iVo,o + Sm^A^i.o + 16m^A/-2,o) 



(197) 



= CACpcfFSa 1 ^^^^ j^^ Zg(TO3,TO3,mD, Wd) + ^ ^ + "lD)Ic(TO3,m3,mD,TOD,0) 



+ [/i(m3) (Ji(mD) + 4rn^/r2(mD)) + h{mT,) {Ji{ms) + 4m^/r2(m3))] 



+ 



(47r)3 6(rf - 1) 



m| m3 + mo ^ rn| m3 + mp _ ^ ^ ^ ' 
ms mD Tod TO3 



(198) 



Cldpgt + 2Cpdp\glg',' + ^g',^ 



-/272 - 27iJi + 2Mo,o 



+ Amlh J2 - 8to^Mi,o + 8m^M2,o] ^^^^ (199) 



= (cYgl + \g'i) dp{dF + 1)A3 [2/2/2 - 4/1 Ji + 8to2/i J2]^^^^ (200) 
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X [Ii{m3)h{m3) - 2Ji(m3) + 4771372(7713)] 



(201) 



^xj i { A = dF(dF + l)'A^/2(m3)/2(m3) (202) 

ij-y-^r^ 1 



i i i i = i;dF{dAh3h{mu) + h'sh{m'a)Yl2{m3) (203) 



) Y = (c;A/l3/l(mD) + /l3-fl(TOD))'^F(dF + l)A32/i(m3)J2(TO3) (204) 



CxdAgi 



-Jlh -hJi+ Mo,o + 2mlhJ2 - 4m^Mi,o + 47n^M2,o 



(205) 



^ X Ci '( >= CA(iA(iF53^3-fi(TO3) [/i(mD)/2(mD) - 2Ji(tod) +4m|,J2(mD)] (206) 

1 '""s/ — 

T X ^ t y A= (o!A/i3/2(mD)+/i3'-r2(mj3))4/2(m3) (207) 

4 V A. ..'V / 
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